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Matrix Elements without Quark Masses on the Lattice 



Oh: 



o 
a^ 
o 

> 

iu: 



X 



A. Donini^, V. Gimenez'^, L. Giusti", G. Martinelli ° 

'^Dep. de Fisica Teorica, Univ. Autonoma de Madrid, Cantoblanco, E-28049 Madrid, Spain. 
^Dep. de Fisica Teorica and IFIC, Univ. de Valencia, E-46100 Burjassot, Valencia, Spain. 
■^Department of Physics, Boston University Boston, MA 02215 USA. 

'^Dip. di Fisica, Univ. di Roma "La Sapienza" and INFN-Sezione di Roma, 1-00185 Roma, Italy. 

We introduce a new parameterization of four-fermion matrix elements which does not involve quark masses and 
thus allows a reduction of systematic uncertainties in physical amplitudes. As a result the apparent quadratic 
dependence of e'/e on ms{fi) is removed. To simplify the matching between lattice and continuum renormal- 
ization schemes, we express our results in terms of Renormalization Group Invariant _B-parameters which are 
renormalization-scheme and scale independent. As an application of our proposal, matrix elements of A7 — 3/2 
and SUSY AF = 2 {F = S,C, B) four-fermion operators have been computed. 



1. Introduction 

Since the original proposals of using lattice 
QCD to study hadronic weak decays sub- 
stantial theoretical and numerical progress has 
been made: the main theoretical aspects of the 
renormalization of composite four-fermion oper- 
ators are well understood the calculation 
of K^-K'^ mixing has reached a level of accuracy 
which is unpaired by any other approach 0|; in- 
creasing precision has been gained in the deter- 
mination of the electro-weak penguin amplitudes 
necessary to the prediction of the CP-violation 
parameter e'/e ptV^'j finally matrix elements of 
AS* — 2 operators which are relevant to study 
FCNC effects in SUSY models have been com- 
puted |^,|o|. Methods, symbols and results re- 
ported in this talk are fully described in |P,|lo|. 

2. Matrix elements vifithout quark masses 

The analysis of K'^ — mixing with the most 
general IS.S = 2 effective Hamiltonian requires 
the knowledge of matrix elements {K^\Oi\K'^) of 
parity conserving parts of the following operators 



01 = ^7^(1-75^^^(1-75)0?'', 

02 = s"(l-75)rf" s^(l-75)d^ 

03 = r(l-75)rf'' s''(l-75)rf", 

04 = r(l-75)rf" s^(l + 75)d^ 



(1) 
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r(l-75)rf'' s''(l+75)d". 



On the lattice, matrix elements of weak four- 
fermion operators are computed from first prin- 
ciples. But, following the common lore, they 
are usually given in terms of the so-called B- 
parameters which measure the deviation of their 
values from those obtained in the Vacuum Sat- 
uration Approximation (VSA). For operators in 
(|^), the i?-parameters are usually defined as 



{K^\0,{^,)\K°) = 
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where _Ci = -5,1,6,2 for {i = 2,..., 5). In 
(||), {K^\Oi\K'^) is parameterized in terms of 
well-known experimental quantities and Bi(fj,) 
{Bk{h) = Biifi)). On the contrary, {K°\0^\K°) 
{i — 2, ... ,5) depend quadratically on the quark 
masses in (^), while they are expected to remain 
finite in the chiral limit and depend only linearly 
on the quark masses. Contrary to fx, Mk, etc., 
quark masses can not be directly measured by ex- 
periments and the present accuracy in their deter- 
mination is still rather poor. Therefore, whereas 
for Oi we introduce Bk as an alias of the ma- 
trix element, by using (^) we replace each of the 
"SUSY" matrix elements with 2 unknown quan- 
tities, i.e. the i3-parameter and mg + rud. To 
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Table 1 

Matrix elements in GeV^ at the renormalization 
scale /i = 2 GeV in the RI scheme obtained with 
the new parameterization and the conventional 



one in 


ref. ^ on 


the same set of data. 






New 


Old 


m 


(3 = 6.0 
this work 


13 = 6.2 
this work 


13 = 6.0 
fl 


P = 6.2 
fl 


(Oi) 
Bi 


0.012(2) 
0.70(15) 


0.011(3) 
0.68(21) 


0.012(2) 
0.70(15) 


0.011(3) 
0.68(21) 


B2 


-0.079(10) 
0.72(9) 


-0.074(8) 
0.67(7) 


-0.073(15)- 
0.66(3) 


-0.073(15) 
0.66(4) 


(03) 
B3 


0.027(2) 
1.21(10) 


0.021(3) 
0.95(15) 


0.025(5) 
1.12(7) 


0.022(5) 
0.98(12) 


(04) 

-04 


0.151(7) 
1.15(5) 


0.133(12) 
1.00(9) 


0.139(28) 
1.05(3) 


0.133(28) 
1.01(6) 


/n \ 


0.039(3) 
0.88(6) 


0.029(5) 
0.66(11) 


0.035(7) 
0.79(6) 


0.029(7) 
0.67(10) 


{O7' ) 


0.019(2) 
0.65(5) 


0.011(3) 
0.38(11) 


0.020(5) 
0.68(7) 


0.014(5) 
0.46(13) 


bT 


0.082(4) 
0.92(5) 


0.068(8) 
0.77(9) 


0.092(19) 
1.04(4) 


0.087(19) 
0.98(8) 


Table 2 
RGI M 
Eq. 


^atrix elements in GeV^ computed as in 
) with as^ ~ . 






/3 = 6.0 


(3 = 6.2 








0.017(3) 


0.016(4) 








-0.051(7) 


-0.048(6) 








0.005(7) 


-0.004(7) 






(0«^^) 


0.072(3) 


0.063(6) 








0.043(3) 


0.032(5) 





overcome these problems, we propose the follow- 
ing new parameterization of AS* = 2 operators 



(i?°|Oi(M)|XO) ^ -Mi/ii?i(^), (3) 

The i3j(/z)-parameters are still dimensionless 
quantities and can be computed on the lattice 
by studying appropriate ratios of three- and two- 
point functions |]lO|. By simply using them, we 
have eliminated any fictitious reference to the 
quark masses, hence reducing the systematic er- 
rors on the corresponding physical amplitudes. 
An alternative parameterization, not used here, 
which can be useful in the future is reported in 



|l0|] . The VSA and B-parameters are also used 
for matrix elements of operators which enter the 
= 1 effective Hamiltonian. Notice that this 
"conventional" parameterization is the only re- 
sponsible for the apparent quadratic dependence 
of e'/e on the quark masses. This introduces a re- 
dundant source of systematic error which can be 
avoided by parameterizing the matrix elements in 
terms of measured experimental quantities and 
therefore a better determination of the strange 
quark mass ms(/i) will not improve our theoreti- 
cal knowledge of e'/e. In this work we have com- 
puted the matrix elements (7r|Of^^|i^) of the four 

fermion operators O^^^ {i — 7,8,9) which con- 
tribute to the A/ — 3/2 sector of e'/e. In the 
chiral limit (tttt can be obtained, using 

soft pion theorems, from {'k~^\0'^/^\K^). For de- 
generate quark masses, = nid = m, and in the 
chiral limit, we find 

lim (^+|0?/V^+) = -M^/.Mim B5(/z) 

m— >0 ' m— »0 

lim (^+|0«^/'|i^+) = -Af,V^ lim S4(^) 

lun {n+\Ol/'\K+) = hl'j^, \im B,{^) . 

3. Renormalization Group Invariant Oper- 
ators 

Physical amplitudes can be written as 

{F\neff\I) = {F\dif,)\I)-C{ti), (4) 

where 0{fi) = {Oi{fi), . . . , Om{ij)) is the operator 
basis (for example the one in (|l|) for the AS' = 2) 
and C(/i) the corresponding Wilson coefficients 
represented as a column vector. C{^) is expressed 
in terms of its counter-part, computed at a large 
scale M, through the renormalization- group evo- 
lution matrix M] 

c{^l) = w[^l,M]c{M), (5) 

where the initial conditions C(Af), are obtained 
by perturbative matching of the full theory to the 
effective one at the scale M where all the heavy 
particles have been removed. Af] can be 

written as (see for example pT| ) 

W[^i, M] = M[n]U[n, M]M-\M] , (6) 
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where U = {as{M)/as{n))^''o'^ ^^'^"^ is the 
leading-order evolution matrix and M(fi) is a 
NLO matrix defined in The Wilson coeffi- 

cients C(/i) and the renormalized operators 0{fi) 
are usually defined in a given scheme, at a fixed 
renormalization scale /i, and they depend on the 
renormalization scheme and scale in such a way 
that only -ffg// is scheme and scale independent. 
To simplify the matching procedure, we propose 
a Renormalization Group Invariant (RGI) defini- 
tion of Wilson coefficients and composite opera- 
tors which generalizes what is usually done for 
Bk and for quark masses. We define 



and using Eqs. (^) and (|^) we obtain 



(7) 



(8) 



The effective Hamiltonian (Q) can be written as 
Heff = O"^^ ■ C^^', where 

C^'Gi ^ ^M]C{M) , 0^^^ = 6i^i) ■ w-^[fi] .(9) 

(jRGi jjy^g^ qRGI ^j^g scheme and scale indepen- 
dent at the order we are working. Therefore the 
effective Hamiltonian is splitted in terms which 
are individually scheme and scale independent. 
This procedure is generalizable to any effective 
weak Hamiltonian. The RGI B-parameters can 
be defined as 



B 



RGI 



(10) 



4. Numerical results 

All details concerning the extraction of matrix 
elements from correlation functions and the com- 
putation of the non-perturbative renormalization 
constants of lattice operators can be found in 
|6lP,p^. In this talk we report the results ob- 
tained in Ip^ . The simulations have been per- 
formed at /3 = 6.0 (460 configurations) and 6.2 
(200 configurations) with the tree-level Clover ac- 
tion, for several values of the quark masses and 
for different meson momenta. The main results 
we have obtained for AS* = 2 and A/ = 3/2 ma- 
trix elements and their comparison with the re- 
sults in |9) are reported in Tables |l| and |[ It is 
interesting to note, as expected from chiral per- 
turbation theory, that matrix elements of AS = 2 



SUSY operators are enhanced respect to the SM 
one by a factor 2 — 12 at /j, = 2 GeV. In Figure |^ 



Figure 1. (^^"1041^^") (GeV) with the new 
(square) and old (circle) parameterization as a func- 
tion of the strange quark mass (in MS at /i = 2 GeV) 
used to obtain the full matrix element from the B- 
parameter. 

we show the strong dependence of {K'^\04\K'^) on 
the strange quark mass when the "conventional" 
parameterization ^ is used, to be compared with 
the result obtained with the new parameteriza- 
tion. The results for the analogous AC = 2 and 
AB = 2 matrix elements presented at the confer- 
ence are reported in iQ . Although we have data 
at two different values of the lattice spacing, the 
statistical errors, and the uncertainties in the ex- 
traction of the matrix elements, are too large to 
enable any extrapolation to the continuum limit. 
For this reason, the best estimate of the central 
values of the B-parameters can be obtained by 
averaging the results obtained at the two values 
of P As far as the errors are concerned we 

take the largest of the two statistical errors. 
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